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A,  INVESTIGATIONS  OF  THE  SOL-GEL  TRANSITION 
A  new  critical  line  in  a  model  for  the  sol-gel  transition 
Site-bond  correlated  percolation  has  been  recently  introduced  by  our  group,  as 
a  model  to  study  solvent  effects  in  the  sol-gel  transition.  This  model  has 
stimulated  further  research  because  of  its  interesting  behavior.  In  particular, 
Delyon,  Souillard  and  Stauffer  have  found  a  line  of  percolation  transitions  without 
the  appearance  of  an  infinite  cluster,  similar  to  the  Griffiths  singularity  in  a 
dilute  ferromagnet.  We  have  complemented  this  study  and  shown  that  there  is  another 
line  of  critical  points  where  the  pair-connectedness  function  becomes  critical 
without  the  disappearance  of  the  infinite  cluster.  The  physical  consequence  of  such 
a  critical  line  is  the  formation  of  a  gel  with  a  network  structure  completely 
different  from  most  physical  gels. 

Migdal-Kadanof f  renormalization  group  analysis  for  a  sol-gel  transition  model 
The  site-bond  Ising  correlated  percolation  model  is  not  easily  analyzed, 
particularly  in  3  dimensions.  This  is  due  to  the  difficulty  in  satisfying  the  Ising 
symmetry  and  the  connectivity  properties.  In  order  to  circumvent  this  problem,  we 
have  carried  out  an  analysis  using  the  Migdal-Kadanof f  R.G.  We  have  proposed  an 
alternative  scheme  which  gives  good  results  for  the  random-site-bond  percolation 
problem.  We  intend  to  generalize  this  procedure  to  the  correlated  case. 

Study  of  droplets  for  a  statistical  mechanical  model  of  gelation 

in  three  dimensions 

Using  a  Monte  Carlo  simulation,  we  have  studied  the  droplet  size  distribution 
for  the  Ising  correlated  site-bond  percolation  problem  In  three  dimensions.  The 
exponent  for  that  distribution  is  found  to  agree  well  with  the  corresponding^'^ 
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exponent  of  the  Ising  model*  We  also  have  discussed  the  relation  between  the 
susceptibility  and  the  second  moment  by  using  the  dilute  Potts  formulation  of  this 
percolation  problem# 


5.  STRUCTURE  OF  BRANCHED  POLYMERS 
Branched  polymer  in  a  solvent 

A  model  for  branched  polymers  in  a  solvent  is  studied  as  a  function  of  three 
parameters:  the  chemical  potential  of  the  sites,  the  chemical  potential  of  the 
bonds,  and  the  nearest-neighbor  interaction  amongst  the  sites#  We  have  shown  that 
this  problem  can  be  mapped  into  a  generalized  percolation  problem.  A  Potts 
formalism  for  this  problem  can  be  formulated  to  which  the  Migdal-Kadanof f  RG  has 
been  applied.  We  have  found  a  surface  of  critical  points  corresponding  to  a 
randomly-branched  polymer  (or  random  animal),  a  line  of  tricritical  points 
corresponding  to  a  polymer  collapse  (0-point),  a  higher  order  critical  point  which 
corresponds  to  percolation,  and  a  first-order  phase  transition  which  corresonds  to  a 
compact  polymer.  The  critical  surface  terminates  at  the  0-line,  and  the  0-line 
ends  at  the  percolation  point.  This  new  prediction  of  a  higher-order  percolation 
point  where  the  branched  polymer  changes  its  structure  is  amenable  to  experimental 
and  computer  verification. 

Flory  theory  for  directed  branched  polymers 
We  have  applied  the  Flory  theory  to  branched  polymeric  systems  in  which  there 
is  a  preferential  direction  for  polymerization  to  occur.  This  anisotropy  gives  rise 
to  new  types  of  critical  behavior  that  can  be  accounted  for  accurately  and  quite 
simply  by  a  Flory  approach.  Specifically,  we  find  an  upper  critical  spatial 
dimension  of  d  *  dc  ■  7  and  correlation  length  exponents  of 
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«  (d  +  11)/ 4 ( d  +2)  9/i'/fd  t 

parallel  and  perpendicular  to  the  anisotropy  respectively. 

C.  NUCLEATIQN  AND  SPINODAL  DECOMPOSITION 

One  of  the  more  interesting  and  important  (from  a  materials  viewpoint)  areas  of 
polymer  physics  is  the  study  of  metastable  polymers.  Both  nucleation  and  growth  and 
spinodal  decomposition  are  important  mechanisms  by  which  polymeric  materials  are 
formed.  We  have  begun  an  integrated  experimental  and  theoretical  program  to 
understand  nucleation  and  growth  and  spinodal  decomposition  in  polymer  systems. 

Spinodals  in  medium-range  systems 

Although  the  spinodal  concept  has  been  questioned,  it  has  been  useful  in 
understanding  experimental  data  in  polymer  systems,  and  we  have  begun  to  study  such 
medium-range  interactions  in  continuum  systems.  No  evidence  of  a  spinodal  was 
found,  but  we  were  able  to  show  that  pseudospinodals  (i.e.,  extrapolations  into  the 
metastable  region)  converge  rapidly  with  the  interaction  range  to  the  true  spinodal 
(if  it  exists).  We  also  observed,  for  the  first  time,  a  breakdown  of  the  classical 
droplet  model. 

The  research  done  so  far  is  in  lattice  systems.  To  make  our  studies  more 
realistic  we  have  begun  to  examine  continuum  systems.  The  interaction  between  the 
particles  is  chosen  to  be  a  hard-core  potential  with  an  additional  attractive  square 
well.  If  the  spinodal  line  exists,  the  compressibility  should  diverge  as  this  line 
is  approached. 

For  short-range  interactions,  as  the  metastable  region  is  probed,  the  system 
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becomes  becomes  unstable  before  the  predicted  spinodal  line  is  reached.  To 
facilitate  a  deeper  probe  into  the  metastable  region,  the  interaction  range  is 
extended  by  increasing  the  width  of  the  attractive  well.  The  compressibili ty  is 
then  measured  as  a  function  of  both  density  and  interaction  range. 

Dynamics  of  supercooled  fluids  and  polymer  melts 
We  have  investigated  the  dynamics  of  supercouled  fluids  and  polymer  melts  based 
on  a  static  mean-field  theory  of  freezing  developed  by  Grewe  and  Klein,  and  the 
dynamical  theory  of  Cahn  and  Hilliard.  A  nonlocal  term  is  introduced  into  the 
standard  Landau-Ginsberg-Wilson  Hamiltonian  and  the  resulting  dynamics  are  explored. 
It  is  shown  that  for  short  times  the  supercooled  liquid  is  unstable  to  density 
fluctuations  of  nonzero  wavevector  k.  The  value  of  k  at  which  this  instability 
first  appears  is  on  the  order  of  the  inverse  range  of  the  interaction  potential. 
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A  statistical- mechanical  model  for  reversible  gelation  is  developed.  This  model  takes 
into  account  solvent  effects,  which  usually  are  neglected  in  the  classical  theory  of  gela¬ 
tion.  The  exact  solution  of  this  model  is  given  for  the  limiting  case  in  which  “loops*’  or 
intermolecular  interactions  may  be  neglected  (Cayley  tree).  The  general  phase  diagram  is 
obtained  and  it  is  shown  that,  with  a  particular  choice  of  a  solvent,  one  can  realize  the  in¬ 
teresting  situation  in  which  gelation  point  and  consolute  point  coincide.  This  point  has 
peculiar  properties  associated  with  the  simultaneous  divergence  of  “connectivity”  and 
thermal  fluctuations.  The  recent  experimental  data  of  Tanaka  and  collaborators  are  in 
good  qualitative  agreement  with  the  predictions  of  the  model. 


L  INTRODUCTION 

Much  of  the  progress  of  the  last  decade  in  sta¬ 
tistical  mechanics  stems  from  the  fact  that  rela¬ 
tively  simple  and  therefore  tractable  models  have 
proved  sufficient  to  describe  extremely  subtle 
cooperative  phenomena.  Three  examples  are 
shown  schematically  in  Fig.  1: 

(i)  A  fluid  near  its  critical  point.  The  Ising  or 
lattice-gas  model1  has  proved  remarkably  success¬ 
ful  in  interpreting  a  wide  range  of  data  near  the 
critical  points  of  fluids.2  This  stems  from  the  fact 
that  the  “essential  physics’*  of  the  problem  is  an 
interparticle  interaction  potential  characterized  by 
a  hard-core  repulsion  and  a  short-range  attraction. 

(ii)  Dilute  and  semidilute  polymer  solutions .  The 
n  =0  limit  of  the  n- vector  model  has  proved  capa¬ 
ble  of  describing  polymer  solutions  in  the  dilute3 
and  semidilute4  regimes,  where  the  magnetic  field 
plays  the  role  of  concentration.5 

( iii )  Polymer  gelation.  The  essential  physical 
feature  of  a  gel  is  connectivity,  and  hence  one  ex¬ 
pects  percolation  phenomena  to  be  relevant.6  As 
we  shall  see,  temperature-dependent  effects  such  as 
those  due  to  the  presence  of  solvent  are  excluded 
from  simple  “pure  percolation."7  It  is  the  purpose 
of  this  paper  to  suitably  generalize  pure  percolation 
in  order  to  incorporate  such  effects.8 

In  Sec.  II  we  shall  describe  our  approach  using 
polyfunctional  condensation,  the  simplest  example 
that  illustrates  the  basic  phenomenon  of  gelation. 

In  Sec.  Ill  we  derive  the  equation  of  state,  while  in 
Sec.  IV  we  describe  the  connectivity  properties. 
Then  in  Sec.  V  we  describe  the  more  general  gela¬ 
tion  of  “vulcanization"  phenomena  in  which  poly¬ 


mers  made  up  of  M  monomers  are  permitted  to 
crosslink.  The  detailed  derivation  of  the  appropri¬ 
ate  formula  are  given  in  Appendixes  A  and  B. 


II.  THE  MODEL  (M  =  1 ) 

We  shall  first  describe  our  model  for  the  sim¬ 
plest  case  of  gelation,  the  polyfunctional  condensa¬ 
tion  of  /-functional  monomers.  Suppose  all  the 
monomers  are  identical,  and  that  each  has  /- 
functional  groups  that  can  react  with  one  of  the  / 
groups  of  another  monomer.  The  simplest  case, 

/  ae  0,  produces  no  reactions  at  all.  The  next  sim¬ 
plest  case,  /« 1,  results  in  dimers  only.  If  /  =  2, 
we  can  have  unbranched  linear  polymers.  For 
/  >  3,  we  form  branched  polymers,  as  illustrated  in 

»  MQOEL 

*  n  *  I 

*  n  *  0 

GEL . *  Gsnn*ci»vity  . *  *  »  l 

[  near  qti  point) 


SYSTEM  - 


ESSENTIAL  physics  - 

Hord  con  fWurtwon  1 
Sncrt-nsnq*  aimxtion  ( 


PLUtO . 

f  nearer i ted  port) 

OtLUTE  POLYMCfl  Ek tried  ***** 


FIG.  1.  Schematic  illustration  of  the  application  of 
specific  model  Hamiltonians  to  capture  the  essential 
physics  embodied  by  various  physical  systems  near  their 
respective  critical  points.  The  symbol  n  refers  to  the 
number  of  components  of  the  order  parameter  in  an  /i- 
vector  model  (isotropically  interacting  n -dimensional 
classical  spins),  while  the  symbol  s  refers  to  the  number 
of  discrete  states  in  a  5-state  Potts  model.  The  “critical 
point**  of  a  dilute  polymer  solution  corresponds  to  the 
limit  iV  —  oo,  where  N  is  the  polymerization  index. 
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(b)  (c) 

FIG.  2.  Illustration  of  the  simplest  gelation  phenomenon,  polyfunctional  condensation  of  /‘functional  monomers. 
The  /-functional  monomer  shown  in  (a)  is  trimethoyl  benzene;  it  has  three  “functional”  groups  which  can  react  to  form 
ether  linkages.  If  /  were  2,  then  the  most  complex  structures  possible  would  be  chains  and  rings.  However,  since  /  >  2 
here,  there  exists  the  possibility  of  forming  branched  networks.  In  (b)  and  (c)  are  shown  beakers  at  successive  stages  of 
reaction.  This  Figure  is  from  Gordon  and  Ross-Murphy  (Ref.  9). 


Fig.  2  for  a  particular  example  with  f  =  3, 
trimethoyl  benzene.9  Each  benzene  ring  has  three 
groups  that  can  react  to  form  an  ester  linkage. 

This  process  is  characterized  by  a  single  parameter 
a,  termed  the  conversion,  which  is  the  fraction  of 
reacted  groups.  Clearly  if  a  =0,  only  monomers 
are  present.  If  0<a  <  1,  there  exists  a  distribution 
of  finite  polymers  of  all  possible  sizes.  However, 
the  probability  of  an  infinite  polymer  or  “gel”  is 
zero  for  all  values  of  less  than  a  critical  value  ac. 
For  a  >  ac,  there  is  a  nonzero  probability  for  the 
occurrence  of  a  single  branched  polymer  that  is  in¬ 
finite  in  spatial  extent.  Thus  the  probability  of  the 
gel  molecule  to  occur  jumps  discontinuously  from 
zero  for  a  <ac  to  unity  for  a  >  ac,  and  hence  the 
connectivity  of  the  system  changes  drastically  at 
a=*ar.  This  ‘‘phase  transition"  is  termed  the  gela¬ 
tion  threshold. 

The  first  successful  model  to  capture  the  essen¬ 
tial  physics  of  the  gelation  threshold  was  proposed 
40  years  ago  by  Flory  and  developed  in  a  series  of 
classic  papers  by  both  Flory10  and  Stockmayer.11 
(AJso  see  the  classic  book  by  Flory.12)  This 
“Flory-Stockmayer'*  (FS)  model  not  only  predicts 
the  occurrence  of  a  gelation  threshold  a=ac,  but 
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FIG.  3.  Phenomenon  of  bond  percolation:  a  Finite 
section  ( 16  X  16)  of  an  infinite  “fence,”  in  which  a  frac¬ 
tion  pB  of  the  links  are  conducting  while  the  remaining 
fraction  =  1  —pa  are  insulating.  Four  choices  of  the 
parameter  pB  are  shown,  (a)  pB  =  0. 2,  (b)  p8  =  0. 4,  (c) 
Pb  =  0. 6,  and  (d)  pB  =  0. 8. 
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FIG.  4.  Analogy  betweeen  the  phase  diagrams  of 
pure  bond  percolation  and  pofyfunctional  condensation 
(gelation  of  /-functional  monomers).  The  role  of  p  in 
gelation  is  played  by  a,  which  is  termed  the  "conver¬ 
sion”  or  "extent  of  reaction”;  it  may  be  thought  of  as 
the  fraction  of  intact  bonds. 


also  provides  precise  predictions  for  the  family  of 
critical-point  exponents  characterizing  the  behavior 
of  various  quantitites  in  the  immediate  vicinity  of 

<V 

For  the  purposes  of  this  paper  it  is  most  ap¬ 
propriate  to  explain  the  FS  model  in  the  context  of 
random-bond  percolation,  which  in  turn  is  illus¬ 
trated  in  Fig.  3.  Suppose  we  have  an  infinitely 
high  and  infinitely  long  wire  fence.  Imagine  also 
that  a  randomly  chosen  fraction  Pb  of  the  links  of 
this  fence  are  conducting  while  the  remaining  frac¬ 
tion  ( 1  —  pB )  are  insulating.  Computer  simulations 
of  a  finite  (16x16)  section  of  this  fence  are  shown 
in  Fig.  3  for  p*=0.2,  0.4,  0.6,  and  0.8. 

Clearly,  for  pB  small,  as  in  Fig.  3(a),  the  system 
consists  of  small  clusters  of  conducting  bonds.  In 
3(b)  the  conducting  fraction  pb  has  doubled,  yet  the 
system  still  consists  of  only  finite  clusters — the 
"scale”  has  increased,  but  not  the  essential  macro¬ 
scopic  conductivity.  In  3(c),  pB  =0.6,  and  the  sys¬ 
tem  is  macroscopically  different:  In  addition  to 
the  finite  clusters,  there  is  a  single  cluster  that  is 
infinite  in  spatial  extent  (of  course,  the  fence  must 
be  infinite  if  the  cluster  is  to  be  infinite).  For 
some  value  of  pB  in  between  Figs.  3(b)  and  3(c), 
there  is  a  threshold  pB\  below  pB  the  fence  cannot 
conduct,  while  above  pB  it  can.  Thus  its  macro¬ 
scopic  properties  change  suddenly  as  a  microscopic 
parameter pB  increases  infinitesimally  from  pB—b 
top$+8  [Fig.  4(a)]. 

Similarly,  below  the  gelation  threshold  ac ,  the 
system  of  Fig.  2  consists  of  only  Finite-size  poly¬ 
mers;  it  cannot,  e.g.,  sustain  a  shear  stress.  Above 
the  threshold  it  can.  Thus  the  macroscopic  prop¬ 
erties  change  suddenly  as  a  microscopic  parameter 
a,  the  extent  of  reaction  (or  equivalently,  the  frac¬ 
tion  of  formed  crosslinks),  increases  infinitesimally 
from  ac  -6  to  at  +6  (Fig.  4(b)], 

The  FS  model  was  formulated  in  a  fashion  that 


at  First  sight  seems  to  be  lattice  independent:  One 
only  requires  that  a  given  polymer  be  forbidden  to 
loop  back  upon  itself.  In  short,  "intramolecular 
interactions"  are  excluded.  This  assumption  is  ful¬ 
ly  equivalent  las  far  as  critical  behavior  is  con¬ 
cerned)  to  the  statement  that  the  polyfunctional 
monomers  be  required  to  occupy  the  sites  of  a 
Cayley-tree  pseudolattice:  To  each  conFiguration 
of  the  /-functional  monomers  there  is  a  one-to-one 
correspondence  with  a  configuration  of  bonds  on 
the  Cayley  tree  with  coordination  number  z  =/. 

The  effect  of  allowing  for  loops  is  clear,  at  least 
on  a  qualitative  level.  Clearly  the  threshold  is  ex¬ 
pected  to  increase,  since  extra  bonds  will  be  formed 
that  will  merely  create  a  loop  rather  than  contri¬ 
buting  to  the  formation  of  an  infinite  branched 
network.  Moreover,  we  expect  that  the  behavior  of 
the  system  in  the  immediate  vicinity  of  the  gel 
point  to  be  characterized  by  different  critical  ex¬ 
ponents.  In  fact,  if  we  are  to  believe  the  utility  of 
lattice  models,  then  it  turns  out  that  exponents  are 
shifted  from  their  Cayley-tree  values  quite  consid¬ 
erably  (Fig.  5).  Of  course,  one  could  well  question 
the  appropriateness  of  a  lattice  model  to  represent 
a  continuum  system.13,14  Hence  much  needed  are 
calculations  for  "continuum  percolation”  that  are 
sufficiently  accurate  to  make  meaningful  predic¬ 
tions  concerning  critical-point  exponents. 

Experimental  evidence  for  departures  for  "classi¬ 
cal"  critical-point  exponents  is  somewhat  incon¬ 
clusive  at  the  present  time.  A  literature  search 
focused  on  this  question  was  recently  earned  out 
by  Brauner,15  who  concluded  that  no  clear-cut 
answer  emerges  despite  rather  extensive  analysis  of 
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FIG.  5.  Illustration  of  the  analogies  between  (a)  an 
ordinary  thermal  phase  transition  (e.g.,  an  Ising  or 
lattice-gas  model),  (b)  bond  percolation,  and  (c)  poly¬ 
functional  condensation. 
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FIG.  6.  Schematic  phase  diagram  of  random-site, 
random-bond  percolation,  where  p,  is  the  site  occupation 
probability  and  pB  is  the  probability  that  a  randomly 
chosen  bond  is  found  to  be  intact. 


existing  data,  A  recent  set  of  experiments  by 
Schmidt  and  Bure  hard  suggests  that  the  Cayley- 
tree  predictions  may  be  quite  adequate,16  Howev¬ 
er,  problems  arise  due  to  the  paucity  of  data  ex¬ 
tremely  close  to  the  gel  point — since  only  there  are 
exponents  expected  to  depan  from  their  classical 
values.17 

In  the  FS  model  solvent  effects  are  not  included. 
Nor  are  temperature  effects  included  in  any  statist¬ 
ical  mechanical  fashion:  All  states  of  a  system 
consisting  of  b  occupied  bonds  are  equally  prob¬ 
able.  This  simplifying  feature  has  great  merit  in 
that  the  FS  model  is  extremely  tractable. 

In  FS  theory  all  sites  are  occupied  by  monomers. 
However,  we  know  that  solvent  effects  are  impor¬ 
tant  in  gelating  systems.  Two  simplifying  assump¬ 
tions  of  the  FS  theory  are  the  following: 

(i)  the  absence  of  solvent  molecules,  and 

lii)  the  absence  of  correlations  between  the  mole¬ 
cules. 

Both  solvent  effects  and  correlations  are  taken  into 
account  in  this  paper. 


A.  Solvent  effects 


responds  to  moving  to  the  right  along  the  horizon¬ 
tal  dashed  line. 

B.  Correlations  among  molecules 

We  shall  assume  that  the  monomers  and  solvent 
molecules  are  not  randomly  distributed  among  the 
sites.  Rather,  we  shall  assume  a  correlation  of  the 
standard  lattice-gas  model  sort  (Fig.  7).  In  specify¬ 
ing  the  interactions,  we  must  consider  that  the 
monomers  can  interact  with  each  other  in  two 
ways.  One  is  the  usual  van  der  Waals  interaction, 
and  the  other  is  a  directional  interaction  that  leads 
to  chemical  bonds. 

The  particle-particle  interaction  of  this  system  is 
reasonably  approximated  by  the  following  nearest- 
neighbor  interactions:  —  the  sol  vent -solvent 

interaction  energy,  —  WaB  the  mono  ^er-solvent 
interaction  energy,  and 


~€bb  = 


-Wbb 

-E. 


(I) 


where  —  WBB  is  the  van  der  Waals  energy  (weight 
pu)  and  —  E  is  the  bonding  energy  (weight  1  —  pu ). 

The  interaction  —€BB  needs  some  further  justifi¬ 
cation.  Two  monomers  can  interact  in  two  dif¬ 
ferent  ways.  The  first  is  the  usual  van  der  Waals 
type  of  attraction,  which  we  have  approximated 
with  a  nearest-neighbor  attraction  —  WBB.  The 


(<j)  (B)  [c; 


Suppose  we  allow  the  sites  to  be  of  two  sorts,  A 
and  B.  A  sites  are  occupied  by  monomers  and  B 
sites  by  solvent.  The  original  “random-bond”  per¬ 
colation  problem  is  now  a  “random-site-bond” 
problem11;  the  FS  critical  point  in  the  simple  phase 
diagram  of  Fig.  4(b)  is  now  an  entire  “line”  of  crit¬ 
ical  points  (Fig.  6).  If  p,  is  the  density  of  A  sites, 
then  FS  theory  corresponds  to  the  special  case 
p,  *  1  (heavy  solid  line).  A  typical  experiment  cor- 


FIG.  7.  (a)  and  (b)  All  sites  are  occupied  by  mono¬ 
mers  (open  circle),  as  in  the  Flory-Stockmayer  model  of 
gelation.  The  wavy  lines  correspond  to  chemical  bonds 
between  two  monomers,  while  p*  is  the  probability  of 
such  a  bond  being  present,  (c)  The  model  of  gelation 
proposed  here.  Each  site  can  be  occupied  by  either  a 
monomer  (circle)  or  a  solvent  molecule  (dot).  The  wavy 
line  corresponds  to  a  reversible  bond.  We  find  />*,  given 
by  Eq,  (12),  is  a  function  of  the  temperature.  The 
monomers  are  correlated. 
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is  the  interaction  which  leads  to  bonds  te.g., 
->jrogen  bonds).  which  we  have  approximated 
A:{h  another  nearest -neighbor  interaction  —  E 
-  WaB .  Of  course  this  second  interaction  oc- 
.un  only  when  the  monomers  are  in  a  particular 
^figuration,  For  a  given  pair  of  nearest -neighbor 
monomers,  the  ratio  between  the  number  of  config¬ 
urations  SE  which  lead  to  hydrogen  bonding  (ener¬ 
gy  -£)  and  the  number  of  configurations  N Ba 
*hich  do  not;  NE/NBB  =(  1  - pu  )/pu  should  there¬ 
fore  be  much  less  than  unity.  Note  that  the  entro¬ 
py  difference  between  the  unbound  state  (with  en¬ 
ergy  -  WBB)  and  unbound  state  (with  energy  —  E) 
is  given  by 

5  =  kBinpu/i\—pu)  . 

The  Hamiltonian  for  a  system  of  SA  molecules 
of  solvent  and  Na  monomers  can  be  written  as 

Uj> 

+  WABintnf+nfnp}. 


(2) 


The  sum  is  over  all  possible  AAf  BB  and  AB 
nearest- neighbor  pairs,  where  fly  =  1  if  site  j  is  oc¬ 
cupied  by  a  solvent  molecule,  and  Uj=0  other¬ 
wise.  Similarly  l\f=  1  if  site  j  is  occupied  by  a 
monomer,  nf =0  otherwise.  Since  each  cell  must 
be  occupied  by  either  a  solvent  or  a  monomer,  we 
have  the  constraint 

n;+nf=i.  oj 

The  partition  function  Zv  j  eBB  |  for  a  given 
configuration  of  interacting  bonds  [  eBB  |  can  be 
written 


2ctpI R.paNa+hbNb 

V,  conf 


-jnNA,NB)]\  ,  (4) 


where  &=  \/kBT,  while  pA  and  pB  are  the  chemi¬ 
cal  potentials  of  species  A  and  B.  The  first  sum  is 
over  all  the  monomer-solvent  configurations  with 
S A  and  Nb  fixed.  The  second  sum  is  over  all 
values  of  N A  —S  —  SB.  By  averaging  over  all 
possible  energies  of  interaction  eBB  one  obtains  the 
partition  function  Zv  for  this  model  of  reversible 
gelation 

Zv  —  Z v ( € BB )  .  (5) 


The  partition  function  (5)  can  also  be  regarded  as 
an  annealed  random-bond  ferromagnetic  problem. 
Using  the  identities 

=  (6) 
j  i 

and  Eq.  (3),  we  can  write  Eq.  (4)  in  the  form 


ZN\eBB\=Z,  2  «P 

m/i 


0 


Mrfr^n  f 

j 


Uj) 


(7) 


Here 


Z,  =exp|  +pAN\  |  (8a) 

is  the  partition  function  of  the  pure  solvent, 

Mrfr =Ms  -M/i  +/(  wab  -  w  aa  >  (8b) 

and 

cl(€bb  )  —  *bb  +  ^ aa  — •  (8c) 

The  sum  in  (7)  is  over  all  possible  configurations 
of  monomers  |  Ilf  (.  Equation  (5)  can  be  written 


2 14  =  Z,  ^  exP 


I  nf\ 


P 


Mefrlnj’+H'X  n?u? 

j  Uj> 


(9) 


where 


e0w=puexp[0a(  WBB )] +( 1  -p*  )exp{0a(£)] 

— P* **!>{#  ^aa  +  —  2 WAB )] + ( 1  —pu  )exp[fftE  +  W m  —  2 WAB )]  .  (10) 


Equation  (9)  is  obtained  using  the  identity 
p.e*"  w»'n‘n‘+{  j  _p> 

(ID 


r - 

In  conclusion,  our  system  is  equivalent  to  a  one- 
component  system  with  an  effective  chemical  po¬ 
tential  Mefr  given  by  (8b)  and  effective  energy  W 
given  by  (10).  From  the  partition  function  one  can 
derive  the  free  energy  and  all  the  thermodynamic 
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properties  of  the  system.  In  particular,  one  can  lo¬ 
cate  the  consolute  temperature  Tc  below  which  the 
system  separates  in  the  different  phases,  and  the 
coexistence  curve. 

One  interesting  question  is  “can  we  derive  the 
gelation  curve  from  the  free  energy?**  That  is,  in  a 
temperature-density  plane,  can  we  derive  from  the 
free  energy  the  curve  which  separates  the  sol  phase 
from  the  gel  phase?  The  answer  is  no  for  the 
model  that  we  are  considering.  In  fact  the  process 
of  gelation  is  related  to  the  connectivity  properties 
of  the  system;  connectivity  properties  have  never 
been  derived  from  a  free  energy. 

A  “gel*’  phase  is  defined  to  be  the  phase  where  a 
nonzero  fraction  of  monomers  are  bonded  together 
via  chemical  bonds  to  form  a  macroscopic  mole¬ 
cule.  In  order  to  calculate  the  gelation  threshold 
6tiT)  we  must  specify  when  a  pair  of  monomers 
are  bonded.  We  require  that  (i)  they  be  nearest 
neighbors  and  (ii)  their  relative  energy  be  —  E. 
Where  two  nearest- neighbor  monomers  satisfy  (ii), 
we  say  that  a  bond  is  present  between  two  mono¬ 
mers.  The  probability  ps  that  such  a  bond  is 
present  between  two  nearest -neighbor  monomers 
can  be  easily  calculated  and  is  simply  given  by19 


Pb  = 


[(l-pM  )e*E] 
[pme*Wu+(l-pu  )e0E] 


(12) 


The  reversible  gelation  that  we  are  describing  is  an 
equilibrium  situation  where  bonds  are  continuously 
created  and  destroyed  on  a  time  scale  short  com¬ 
pared  to  most  observations  (weak  gels).  It  definite¬ 
ly  does  not  apply  to  strong  gels.  Thus  the  time 
scales  of  weak  gels  are  analogous  to  those  of  an¬ 
nealed  random  magnets,  while  strong  gels  are 
analogous  to  quenched  random  magnets.20 

The  problem  of  calculating  <$f(D  is  in  some 
respects  analogous  to  the  usual  site-  or  bond- 
percolation  problem.7  However,  it  is  more  com¬ 
plex  for  the  following  two  reasons:  (a)  In  the 
“pure**  site-percolation  problem,  the  particles  are 
randomly  distributed,  while  here  they  are  correlat¬ 
ed  according  to  the  Hamiltonian  (2)  (“correlated” 
percolation),  (b)  In  site  percolation,  the  vertices 
can  be  occupied  or  not  and  the  bonds  are  always 
present.  In  bond  percolation  the  vertices  are  all 
occupied  by  the  particles  and  the  bonds  may  be 
present  or  absent.  In  our  model  the  vertices  may 
be  occupied  or  not  and  also  the  bonds  may  be 
present  or  absent  (“site-bond”  percolation).  The 
correlated  percolation  problem21-25  and  the  site- 
bond  percolation  problem18  have  each  been  treated 
separately,  while  we  have  proposed  treating  both 


problems  simultaneously.4 

From  the  above  considerations,  n  follows  that  if 
Pb  < Pc*  where  pc  is  the  pure  bond- percolation 
threshold,  there  is  no  gelation  no  matter  how  high 
the  monomer  density.  Therefore  from  Eq.  (12) 
there  exists  a  limiting  value  of  the  temperature 
Tm^>  with  Pb(T< nu)  ~pc  above  which  there  is  no 
gelation.  From  Eq.  (12)  it  follows  that  Tm„  does 
not  depend  on  the  nature  of  the  solvent. 

Thus  far  we  have  presented  a  model  to  describe 
the  sol-gel  phase  transition  for  weak  gels.  An  ex¬ 
plicit  expression  for  pB  as  function  of  temperature 
has  been  obtained  and  a  value  of  temperature  T mA1, 
independent  of  solvent,  is  predicted  above  which 
there  is  no  gel.  In  the  next  two  sections  we  will 
solve  this  model  for  the  interior  of  the  Cayley  tree. 
We  will  first  calculate  the  equation  of  state  (Sec. 
Ill)  and  then  the  connectivity  properties  (Sec.  IV). 
Results  on  this  lattice  correspond  to  making 
Flory’s  assumption  of  no  intramolecular  interac¬ 
tions  and  hence  the  closed-form  expressions  we  ob¬ 
tain  provide  a  useful  anchor  point  for  theoretical 
descriptions  of  this  model  system. 


III.  EQUATION  OF  STATE  (A#  « l ) 

In  the  preceding  section  we  have  shown  that  the 
monomer-solvent  solution  is  equivalent  to  an  effec¬ 
tive  one-component  system  with  the  partition  func- 


(o ) 


FIG.  8.  (a)  Cayley  tree  of  coordination  number  /  ^  3. 
The  center  0  is  connected  to  branches  with  on  gin  in 
1,2,3.  (b)  Example  of  branch  with  origin  in  0. 
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non  given  by  Eq.  (9).  For  such  a  system  we  calcu¬ 
late  here  the  equation  of  state  for  the  interior  of 
the  Cavley  tree,  following  closely  the  procedure  of 
Ref  24. 

To  fix  the  ideas,  consider  the  Cayley  tree  or 
Bethe  lattice  of  coordination  number  (“functionali¬ 
ty*!  /  =  3,  with  a  “center'’  denoted  0  [Fig.  8(a)]. 

Any  site  is  connected  to  /  branches.  Any  bond,  if 
jut  divides  the  tree  in  two  branches.  For  example, 
if  we  cut  the  bond  [03]  we  have  one  branch  with 
origin  0  [Fig.  8(b)]  and  another  with  origin  at  3. 

We  split  the  partition  function  Zy  *Z.V(  4- ) 

—  1  where  Zy(  4- )  [Zy(  —  )]  is  the  partition 
function  under  the  condition  that  site  0  is  occupied 
empty).  They  satisfy  the  following  relations: 

Zy (  4*  )» ZyQ(  4*  )[e*WZNy{  4-  )  +  Zy3(  —  )]  , 

(13a) 

Z y (  -  )=ZyQ(  -  )[Zy3(  4*  >  +  Z*3(  -  )]  ,  (13b) 

where  Zy  (  4- )  [ZyQ(  —  )]  is  the  partition  function 
of  the  branch  0  [Fig.  8(b)]  under  the  condition  that 
site  0  is  occupied  (empty).  Analogous  definitions 
hold  for  the  partition  function  relative  to  branch  3. 
.V0  and  jV3  are  the  number  of  sites,  respectively,  in 
branch  0  and  3.  The  density  of  monomers  ^  is 
clearly 

4  =  lim  [Zy(  4"  )]/[Zy(  +  )  +  Zy(  —  )]  •  (14) 

.V— •• 

From  Eqs.  (13)  and  (14) 

<i=b(V*'+i)l/b’2^','+2j'  +  i] .  (15) 

where 

y=  Hm  [Z*0(  +  )]/[Z*0(-» 

•Vq  ® 

-  lim  [Z„A  +  )]/[Z„yi-)].  (16) 

.v3— « 

Here  we  have  made  the  assumption  of  translational 
invariance.  The  translational  invariance  condition 
is  equivalent  to  neglecting  the  surface  of  the  Cay¬ 
ley  tree. 

We  need  an  equation  for  y.  To  do  so  we  write 
the  following  relations: 

2 y  ( 4* )»***•*  17  [e*WZN{ 4- )+Zy  (  — )]  t 

(l  7a) 

zv  -  >-  + >+**,<  - )] .  d7b) 

tm  1 

Here  we  have  introduced  the  partition  functions 


for  the  branches  1  and  2  which  have  been  obtained 
by  cutting,  respectively,  bonds  [01]  and  [02]. 

Dividing  Eq.  (17a)  by  Eq.  (17b)  and  using  the 
translational  invariance  assumption,  we  obtain  in 
the  thermodynamic  limit  the  desired  expression  for 

y> 

y  ):Ay+  I  )*  .  1 18) 


A.  General  functionality 


First  we  introduce  the  new  variables 


z—e~ew/1 , 

1 19a) 

II 

a 

9 

£ 

1 

(19b) 

P\=y/z  . 

(19c) 

Equations  (15)  and  (18),  generalized  to  arbitrary 
functionality  /,  become 

^  =  (|ii+/i,z)/(fi?  +  2/i,r  + 1)  , 

(20a) 

p=P\[(p\Z  +  \)npx+z)Y~'  , 

(20b) 

which  are  the  desired  equations  of  state  coinciding 
with  the  Bethe  approximation.24  From  (19a),  (19b), 
and  (8b),  /i  can  be  written  in  the  following  way: 

fi=cxp(0\  [MB-FMs(0)]-[/iw+^(0)J  ))  , 

(21a) 

where (0)  =  {  \)fWBB  and,u*(0)  =(  \)P#aa  arc 
the  chemical  potential  of  species  B  and  A,  respec¬ 
tively,  in  the  absence  of  species  A  and  B.26  Here 

e'*"  =pueBWl,  +  {  1  -pu  )e»E  .  (21b) 

From  (20a)  it  follows  that  p  =  1  corresponds  to 
the  disordered  phase.  The  consolute  temperature 
Tc  (or  equivalently  zc )  and  the  consolute  density  <£f 
are  obtained  from  the  equations: 


FIG.  9.  (a)  Coexistence  curve  for  the  monomer  case 
(M  =  1)  and  (b)  for  the  polymer  case  [M  >  >  1). 
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FIG.  10.  Example  of  ghost  sites  for  a  section  of  the 
square  lattice  (a)  and  Cayley  tree  ib). 


extended  by  Conigho22,24  for  correlated  percolation. 
Here  we  calculate  P  in  the  presence  of  the  "ghost” 
field  h  (Refs.  28  and  29)  which  plays  the  same  role 
as  the  magnetic  field  in  a  ferromagnet.  To  be 
more  precise/9  we  introduce  a  ghost  site  which  is 
connected  to  every  site  with  probability  h  (see  Fig. 
10).  Consequently  all  the  sites  that  are  connected 
with  the  ghost  sites  are  connected  to  an  infinite 
cluster.  The  advantage  of  introducing  the  ghost 
field  stems  from  the  relation2* 


or 


(22a) 


dP 

a*  Lo' 


(27) 


and 

=0  or  “-^-=0.  (22b) 

o*  dMl 

From  Eqs.  (20a),  (20b),  (22a),  and  (22b)  we  have 

*(/— 2)//,  <*,  =  7  .  (23) 

The  spinodal  curve  <bt  is  obtained  from  Eqs. 

(20a)  and  (22a), 


Starting  from  the  origin  of  the  elementary  cell 
[Fig.  8(a)]  there  are  /  branches  that  emanate  from 
the  origin.  Given  that  the  origin  is  occupied  by  a 
monomer,  let  Q  be  the  probability  that  moving 
along  one  of  the  branches  there  is  no  infinite  clus¬ 
ter  attached  to  the  origin.  The  probability  P  that 
the  origin,  occupied  by  a  monomer,  belongs  to  an 
infinite  cluster  is  given  by 

/>«1_U-A){?>\  (28) 


+  +  •  (24) 

Here  is  the  solution  of  Eq.  (22a)  and  is  given 
by 

Mu=4f±(7F2- l)l/2  (23a) 

and 

F  =*-'</ -2-/zJ)  .  (25b) 

The  +  and  —  signs  in  (23a)  correspond  to  the 
two  different  branches  of  the  spinodal  curve.  The 
coexistence  curve  is  obtained  by  putting  p  =  1  and 
solving  (20a)  and  (20b)  for  as  function  of  T 
(Fig.  9). 

IV.  CONNECTIVITY  PROPERTIES 

In  this  section  we  investigate  the  connectivity 
properties  of  the  system.  In  particular,  we  are  in¬ 
terested  in  the  following  quantities:  (i)  the  percola¬ 
tion  probability  F,  defined  to  be  the  probability 
that  a  given  monomer  belongs  to  the  infinite  clus¬ 
ter,  (ii)  the  mean  cluster  size  5, 

•  (26) 

Here  n,  is  the  average  number  of  clusters  (mole¬ 
cules)  per  site  of  s  monomers. 

To  calculate  P  we  follow  closely  the  formalism 
introduced  by  Essamr  for  random  percolation  and 


Clearly  Q  satisfies  the  following  relation 

Q  =  \-apB+apB(  \  ,  (29) 

where  a  is  the  probability  that  one  of  the  peri¬ 
pheral  sites  of  the  elementary  cell  [e.g.,  site  1  of 
Fig.  8(a)]  is  occupied  under  the  condition  that  the 
origin  is  occupied, 

a  =  <n£nf>/<n?> ,  oot 

where  the  angle  brackets  stand  for  the  usual  sta¬ 
tistical  average 

<..•>  =  2  01) 

co nf  con^ 

Now  a  has  already  been  calculated,24  with  the  re¬ 
sult 

a  =}i\/(p\+z)  ,  (32) 

where  l  is  given  by  Eq.  (19a)  and  by  Eqs.  (1 9c) 
and  (20b)  (see  also  Appendix  A  for  a  more  general 
case  1).  pB  is  the  probability  that  the  two 
monomers  at  the  origin  0  and  at  the  peripheral  site 
1  are  in  a  “bound  state,”  and  is  given  by  Eq.  (12). 
From  (27)  — (29),  we  have 

0  (\+apMQf~~l) 

S~{\-(f-\)aptQ'-'')  '  (33' 

where  Q  is  calculated  for  h  =0. 
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FIG.  1 1.  Coexistence  curve  for  the  monomer-solvent  binary  mixture  and  sol-gel  phase  boundary  for  three  different 
solvents  (M  =  1).  The  solvents  have  been  chosen  in  such  a  way  that  the  consolute  point  is  (a)  in  the  sol  phase  (b)  on  the 
gelation  curve,  and  (c)  in  the  gel  phase.  Tp  is  the  temperature  at  which  the  gelation  curve  meets  the  coexistence  curve. 
The  situation  depicted  in  Ic)  is  not  expected  to  be  observable  for  usual  experimental  time  scales. 


We  note  that  (29)  is  the  same  equation  as  for 
random  site  percolation.  In  this  case  apB  becomes 
the  monomer  density  The  equation  for  h  =0 
has  been  studied  previously.27 

The  solution  near  the  percolation  threshold  is 


Qix)=  1,  x  <xc  (34a) 

CU)~l-[2/(/-2)]e,  x>*f.  (34b) 

Here  xe  =  (/-l)",1  x  =apBt  and  e  =  U  -xe)/xe. 
From  Eqs.  (28),  (29),  and  (33),  it  follows  that  for 

X  >XCy 

P~[2f/(f-2)]€  (35a) 

and 

Sa// 1*|  .  (35b) 

The  equation  for  the  gelation  threshold  <fig  is  deter¬ 

mined  by  the  condition 

=(/  - 1 P1  •  (36) 

From  Eqs.  (20a)  and  (32)  we  find 

+r2[2(/  —  l)p*  —  1] )  ,  (37) 

which  is  identical  to  the  critical  threshold  in  the 


pure  correlated  case.  The  only  difference  is  the 
change  of  (/— !)-*</  — Dp*.  We  see  that  the 
bond  dilution  reduces  the  effective  connectivity  of 
the  lattice.  We  also  note  that,  consistent  with  our 
general  result,  we  find,  from  (36)  and  from  the  fact 
that  z  <  1,  the  limiting  value  of  the  temperature 
T m4*  above  which  there  is  no  gelation: 

PidW-tz-ip1 .  (38) 

Note  that  (/-l)"1  is  the  percolation  threshold  for 
the  pure  random  bond  percolation. 

In  Fig.  1 1  we  have  plotted  the  sol-gel  phase 


boundary,  Eq.  (37),  together  with  the  coexistence 
curve  for  the  binary  mixture  of  monomers  and  sol¬ 
vent,  for  three  good  solvents.  Changing  the  sol¬ 
vent  corresponds  to  changing  the  parameters 
WBB  —  2Wab.  The  solvent  parameters  have  been 
chosen  in  such  a  way  that  the  consolute  point  is  in 
the  sol  region  [Fig.  11(a)],  on  the  gelation  curve 
[Fig.  11(b)],  and  in  the  gel  region  [Fig.  11(c)].  We 
stress  two  interesting  features: 

(i)  For  all  solvents  there  is  a  temperature  Tp 
(below  the  consolute  temperature  Tc )  at  which  the 
coexistence  curve  crosses  the  gelation  curve.  For 
T  <  Tp  we  have  coexistence  between  sol  phase  and 
gel  phase.  In  addition,  in  Fig.  1 1(c)  for  Tp<T 

<  Tc  we  have  two  possible  gel  phases.  Note  that 
in  a  real  system,  this  situation  is  very  difficult  to 
realize  because  the  monomers  form  very  large  clus¬ 
ters  for  short  time  scales  when  the  demixion  con¬ 
centration  <£f  is  near  the  gel  curve.  Therefore 
tends  to  decrease.  Hence  the  situation  depicted  in 
Fig.  1 1(c)  is  nos  expected  to  be  realized  in  prac¬ 
tice.30  Only  for  an  infinite  time  scale,  the  mono¬ 
mers  act  as  single  elements  and  phase  separation 
would  occur  at  ^c. 

(ii)  By  changing  the  solvent  properties  it  is  possi¬ 
ble  to  realize  the  interesting  case  in  which  the  con¬ 
solute  point  lies  on  gel-sol  phase  boundary,  as 
shown  in  Fig.  11(b).  We  find  that  this  special 
point  Q  is  realized,  if  for  fixed  values  of  the  solute 
parameters  (pu,EyWBB)  the  solvent  parameters 

( ~  2  Wab  )  and  Tc  are  related  by 

-  wnkTe  ,  . 

pB  =  \—e  f  3 1  —  zc  .  (39a) 

In  fact,  from  (37)  using  (39a)  and  z  =zc  =// 

(/  -2)  we  find  *  y.  This  means  that  critical 
point  and  gelation  threshold  coincide.  Equation 
(39a)  is  actually  valid  for  any  lattice  and  for  ail 
d.31 

This  particular  point  Q  is  characterized  by  the 
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divergence  of  two  lengths.  One  is  the  usual  corre¬ 
lation  length  §  which  diverges  at  the  consolute 
point,  and  the  other  is  the  characteristic  linear  di¬ 
mension  of  the  finite  clusters  (the  “connectedness 
length*')  which  diverges  on  the  sol-gel  phase 
boundary.  It  is  interesting  to  study  the  nature  of 
point  Qy  which  in  some  respects  is  analogous  to  the 
point  (7*  ~0,p  ~pc )  in  the  T  —p  phase  diagram  of 
a  randomly  dilute  ferromagnet. 

From  Eqs.  (35a)  and  (35b)  we  see  that  along  all 
the  gelation  curve,  including  the  Q  point,  the  per¬ 
colation  critical  exponents  are  the  classical  ones  if 
the  parameter  x  =apB  is  chosen  as  variable.  Ex¬ 
perimentalists,  however,  do  not  measure  exponents 
along  a  path  of  varying  x  but,  for  example,  along 
the  T  axis.  Therefore  let  us  choose  the  value  of 
Pb*  given  by  Eq.  (39a),  so  that  we  obtain  the  situa¬ 
tion  in  which  the  Q  point  occurs.  We  can  ap¬ 
proach  the  Q  point  along  the  coexistence  curve 
ifj.  =  1 )  using  as  variable  |  T  —  Tc  | .  We  find  criti¬ 
cal  exponents  different  from  those  that  we  find  if 
we  approach  the  Q  point  from  T  >TC.  In  order  to 
show  this  we  must  calculate  near  Te 

x  —xc=pBa(fiui)—(f  —l)~]  ,  (39b) 

with  pB  given  by  Eq.  (39a),  a{puz)  given  by  Eq. 
(32),  and  /i,  from  Eq.  (20b)  with  /is  l .  It  is 
straightforward  but  tedious  to  show  that  x  —  xe 
^T~Te  for  T  >TC  and  x  — —  Dl  72  for 
T  <Te  along  the  coexistence  curve.  Consequently 
the  mean  cluster  size  diverges  with  a  classical  ex¬ 
ponent  yp  =  1  for  T  >  Tc  and  with  a  different  ex¬ 
ponent  yp  =  y  for  T  <  TV . 

This  asymmetry  above  and  below  Te  is  not 
found  in  low  dimensionality  using  a  Migdal- 
Kadanoff  renormalization-group  approach.31 
However  the  critical  exponents  are  found  to  be 
given  by  random  percolation  exponents  along  the 
gdBtion  curve,  with  a  crossover  to  different 
behavior  as  the  Q  point  is  approached.  This  same 
kind  of  crossover  has  also  been  found  near  six  di¬ 
mensions  using  the  ^-expansion  technique,32  where 
the  calculations  have  been  done  only  in  the  disor¬ 
dered  phase  ( T  >  Tt;  p,  =  1 ).  More  study  therefore 
is  required  to  fully  understand  the  critical  behavior 
near  the  Q  point  especially  as  a  function  of  the 
dimensionality  of  the  system. 

V.  POLYMERS  IN  A  SOLVENT  [M  >  \ ) 

In  this  section  we  want  to  extend  the  previous 
treatment  to  the  case  of  a  solution  of  polymers  in  a 
solvent.  Treatment  of  chain  percolation  can  also 


FIG.  12.  Polymer  chains  embedded  on  a  square  lat¬ 
tice.  The  wavy  line  is  a  bond  between  two  monomers  of 
two  different  chains. 


be  found  in  Refs.  5  and  33.  The  polymer  is  made 
of  A/  monomers,  each  one  of  which  occupies  a  lat¬ 
tice  site  (Fig.  12).  The  solvent  molecules,  as  in  the 
case  Msl,  fill  the  remaining  sites  of  the  lattice. 
The  monomer-monomer,  solvent-solvent,  and 
monomer-solvent  interactions  are  given  by  Eq.  (1) 
as  before.  The  partition  function  becomes 
Zjv  jc**  I  for  a  given  configuration  of  monomer  in¬ 
teraction  |6j£  (  can  be  written  as: 

Zj V  I  €BB  !  t  *BB  !  ♦  <40a) 

*  »0 

Here  Z,  is  the  partition  function  of  the  pure  sol¬ 
vent,  given  by  Eq.  (8a),  while  Zn  is  the  contribu¬ 
tion  to  the  partition  function  which  comes  from 
those  configurations  with  n  polymers, 

-«*■"*  2  (40b) 

config 

The  sum  is  over  all  configurations  of  the  n  poly¬ 
mers,  SBB  is  the  number  of  nearest -neighbor  pairs 
of  monomers,  and 

MdT=  Mpoly- V/i, 

—  [(/  — 2)Af+2](  WM-WAB) 

-(M-\)WAA.  (41) 

where  and  are  the  chemical  potentials, 
respectively,  of  the  polymer  and  the  solvent.  The 
origin  of  (41)  can  be  easily  understood  considering 
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FIG.  13.  Examples  of  configurations  of  two  dimers 
without  nearest-neighbor  monomers  and  lb)  with  two 
nearest- neighbor  monomers. 

the  example  of  two  dimers  in  =  2 M  =2),  as  illus¬ 
trated  in  Fig.  13.  There  are  two  possible  configu¬ 
rations.  We  have  to  calculate  the  energy  and  the 
chemical  potential  of  a  system  made  of  ( iV  —  4 ) 
solvent  molecules.  Since  the  contribution  to  the 
chemical  potential  is  nAN  in  Z„  we  must  subtract 
such  a  term  and  add  n^poiy  due  t0  l^e  presence  of 
two  polymer  molecules.  The  contribution  to  (41) 
due  to  the  chemical  potential  is 

.  {42) 

The  contribution  to  (41)  due  to  the  interaction  is 

Wil'mnUf-2)M+2]lWAA-WMB) 

(43) 

The  sum  l^111  =nfieff  +  proves 

Eqs.  (40)  and  (41). 

Taking  in  (40b)  the  average  over  all  possible 
values  of  €bb*  the  partition  function  Z^=Z^  j  €bb  I 
becomes 

(44a) 


Zfi  . 


Here 


conf 


(44b) 


where  W  is  given  by  Eq.  (8c). 

In  conclusion,  the  partition  function  of  our  sys¬ 


tem  is  equivalent  to  the  grand  partition  function  of 
a  system  of  polymers  with  an  effective  chemical 
potential  given  by  Eq.  (41)  and  an  effective 
monomer-monomer  interaction  W  given  by  Eq. 

(10). 

As  for  the  case  M  =  1,  we  derive  the  equation  of 
state  in  closed  form  for  the  intenor  of  the  Cayley 
tree.  The  equation  of  state  is  derived  in  Appendix 
A.  Here  we  give  only  the  result  for  the  monomer 
density 


<b  — R  [fi]  +  fi\Z) /[Rfi]  +  {\  +  R)fX\Z  4-1] 


Here 


R  = 


JiLzlL 


[(/-!)(/  — 2) +  2</-l)/Af]  ’ 


and 


= exp  |  01/ipoiy  +Mpolyl 0  >  1 

-M[fiA  +^(0)11 


(45) 

(46a) 

(46b) 

(46c) 


where  ^^(0)  is  the  chemical  potential  of  the  po¬ 
lymers  in  absence  of  solvent,  and  ^(0)  is  the 
chemical  potential  of  the  solvent  in  absence  of  po¬ 
lymers;  explicit  expressions  are  in  Appendix  A  and 
z  has  been  defined  in  Eqs.  (19a)  and  (10).  If  we  set 
ilfsl,  Eqs.  (45)  and  (46)  reduce  to  Eqs.  (20a)  and 
(20b)  derived  in  Sec.  III. 

The  consolute  temperature  Tc ,  or  equivalently 
zcy  is  obtained  from  the  two  equations 


-^=0  (or  equivalently,  =0)  , 
—^7  =0  (or  equivalently,  =0)  . 

d<*2  atf 

From  Eqs.  (45)  and  (47)  we  find 


i/2 


(47a) 

(47b) 

(48a) 


and 

4c_(£  +RM'/2ze)/[R  +(l+R)MW2zc+M]  . 

(48b) 

Clearly  for  large  Af, 

re-(/-I)*l[(/-l)(/^2)l,/2 
and  d>e  -  RzcM~w2. 


The  spi nodal  curve  <t>s  can  be  obtained  from  Eqs. 
(45)  and  (47a) 

<bs  ~R  i^L\s  +ii\sz)/[Rfi}s  4-(  1  +R)fi\sz  4- 1  ]  ,  (49) 

where  fi]s  is  the  solution  of  Eq.  (47a)  and  is  given 
by 

Mi,  =  ^/2±(F2/4-Af',),/2 


(50a) 
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—  —  \  )  +  l]z2-A/(/  -2)\/Mz 

(50b) 

For  the  case  M  =  l,  in  order  to  Find  the  coex¬ 
istence  curve  d>cocx,  it  was  enough  to  put  ji  =  1  be¬ 
cause  of  the  symmetry.  This  time  we  have  to  ap¬ 
ply  the  “equal  area”  rule.  The  conditions  for  the 
coexistence  curve  are  that  the  pressure  P  and  the 
chemical  potential  must  be  the  same  in  coexisting 
phases  A  and  5.  From  the  relation  dP/d^cff=p, 
where  p  is  the  density  of  chains  related  to  the  den¬ 
sity  of  monomers  by  <6=Afp.  The  condition  for 
the  pressure  to  be  the  same  is 
rB  rB 

JA  &  =  JA  P^Mefr=0 

or 

<£d(lnjx)  =0  , 

and  for  the  chemical  potential  p(/0=ju(P).  In 
Fig.  14  we  have  plotted  schematically  as  a 
function  of  T.  Note  that  the  critical  density  <t>c  is 
located  much  nearer  to  zero  than  in  the  symmetric 
case  M  =  1  (Fig.  9). 

The  connectivity  properties  of  the  system  of 
chains  can  also  be  derived.  Here  we  give  only  the 
result  for  the  gelation  curve  <£f  (the  derivation  is 
given  in  Appendix  B), 

*f=*z2(/-l)/t  [U +/*)</- 1) 

—  l]r2-M/  — 2)2 )  .  (51) 

Here 

f-\=z[M{f-\)  +  \-M]pB  ,  (52) 

and  Pa  is  given  by  Eq.  (20b).  The  discussion  of 
the  results  for  St  >  1  is  identical  to  the  M  =*  1  case, 
and  the  results  are  indicated  in  Fig.  14. 


In  conclusion  we  have  proposed  a  model  which 
applies  to  weak  gels.  We  have  solved  the  model  in 
closed  form  for  the  Cayley  tree.  The  general  phase 
diagrams  are  in  qualitative  agreement  with  the  ex¬ 
perimental  data  of  Tanaka  et  a/.54  and  Ruiz 
et  ai 35  This  model  can  also  be  adapted  to  explain 
peculiar  effects  inherent  in  the  At  >>  1  system.36 
An  interesting  situation  occurs  when  for  a  particu¬ 
lar  solvent  the  gelation  curve  ends  at  the  consolute 
point.  This  point  is  a  higher-order  critical  point 
where  both  the  correlation  length  and  the  connect¬ 
edness  length  diverge. 

In  the  Cayley-tree  solution,  which  is  valid  for 
high  dimensionality,  the  critical  behavior  in  the 
sol-gel  transition  exhibits  an  asymmetry  above  and 
below  the  consolute  temperature.  This  asymmetry 
is  not  found  in  a  renormalization-group  approach 
in  low  dimensionality.31  More  study  needs  to  be 
done  to  investigate  the  nature  of  this  Q  point  as 
function  of  the  dimensionality.  This  point  presents 
analogies  with  the  end  of  the  line  of  critical  points 
in  the  dilute  ferromagnets.  Such  lines  of  critical 
points  cannot  be  obtained  with  the  previous  simpli¬ 
fied  theory  of  gelation. 

After  this  work  was  completed,  several  exten¬ 
sions  of  the  present  model  were  developed.  In  par¬ 
ticular,  Barrett37  has  recently  generalized  the 
present  model  to  incorporate  features  displayed  by 
the  polymeric  system  hyaluronic  acid,  which  un¬ 
dergoes  an  order-disorder  transition.  There  are 
two  types  of  side  groups,  amide  and  carboxyl,  as 
well  as  two  types  of  “solvent  molecules,’*  potassi¬ 
um  ions  and  phosphate  ions.  Very  recently,  De- 
lyon  et  a/.38  found  that  the  present  model  displays 
a  peculiar  transition  in  the  region  below  Tf .  A 
physical  interpretation  of  this  effect  has  been  sug¬ 
gested  by  Klein  and  Stauffer.39 
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FIG.  14.  Same  as  Fig.  1 1  except  that  M  is  large  (polymer). 
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FIG.  15.  Polymer  chains  (heavy  line)  embedded  on  a 
Bethe  lattice  of  coordination  number  /  =  3  The  open 
circles  are  the  monomers.  The  dots  are  the  solvent  mol¬ 
ecules.  The  wavy  lines  are  the  bonds  between  monomers 
of  different  chains. 

The  special  case  of  the  present  model  in  which 
bonds  are  considered  broken  with  probability 

q± b  ~  1  —Pb  =  exp(  —U/kT) 

was  introduced  by  Coniglio  and  Klein31  as  a  candi¬ 
date  for  representing  “droplets”  at  the  critical 
point  of  an  Ising  model  with  coupling  constant  J. 
Numerical  studies  in  two  and  three  dimensions  in¬ 
dicate  that  these  site-bond  percolation  clusters  in 
fact  do  behave  just  like  Ising  droplets  in  that  their 
spatial  extent  becomes  infinite  at  the  critical 
point.40  Thus  we  see  that  a  model  introduced  in 
connection  with  polymer  gelation  would  appear  to 
be  of  relevance  for  a  completely  different 
problem — that  of  defining  Ising  droplets  near  the 
critical  point  of  a  thermal  phase  transition. 
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APPENDIX  A:  POLYMERS  IN  A  SOLVENT, 
EQUATION  OF  STATE 

Here  we  derive  the  equation  of  state  of  a  solu¬ 
tion  of  polymers  in  a  solvent  for  the  interior  of  the 
Cayley  tree.  The  result  was  given  in  Sec.  V,  Eqs. 
(45)  and  (46).  The  polymer  is  made  of  Xf  mono¬ 
mers,  each  one  of  which  occupies  a  lattice  site 
(Fig.  15).  The  partition  function  Z„  is  given  by 
Eqs.  (44a)  and  (44b)  which  we  recall  here  for  con¬ 
venience 

zv=z,  2Z»*  (Al) 

n  -0 

with 

.A2> 

cxmfiz 

All  the  other  quantities  are  defined  in  Sec.  V. 

Consider  first  the  case  M  —  2.  Following  the 
derivation  for  the  M  =  1,  we  consider  the  Cayley 
tree  of  coordination  number  3  [Fig.  8(a)]  and  wnte 

Zs  =  Z#{  — )  q- Zy(  4- )  »  (A3) 

where  Z^(  — )  and  Zv(  4- )  are  the  partition  func¬ 
tions  under  the  condition  that  the  origin  0  is, 
respectively,  empty  or  occupied  by  a  monomer.  In 
the  latter  case  we  can  write 

Zy (  4* )  =  Zv[0, 1  ]  4-  Zv[0,2]  4-Z\[0, 3]  ,  (A4) 

where  Z.v[0,/]  U  =  1,2,3)  is  the  partition  function 
under  the  condition  that  the  dimer  occupies  bond 
[0,/].  The  following  relations  hold; 

Z>[0,l]=e*‘rfrri(e,,'Zivi<  +  )Zsi  -  I]  , 

iM2  (A5) 

3 

£*(  —  )=  f7[Zv((  4“  )4-Zv  (  —  )]  •  (A6) 

l-l 

Here  Zv,(  —  )  and  ZsS  4- )^ZV,[2,6]  +  ZVJ2,7] 
are  the  partition  functions  of  branch  2  under  the 
condition  that  site  2  is,  respectively,  vacant  or  oc¬ 
cupied  by  a  monomer.  Zv,[2,6]  and  Zv,[2,7]  are 
the  contributions  relative  to  the  configurations  in 
which  the  dimer  occupies,  respectively,  the  bond 
[2,6]  and  [2,7].  Analogous  definitions  hold  for  any 
other  Z.v  (  4- )  and  Zv  (  —  )  and  analogous  relations 
hold  for  Zv[0,2]  and  Zv[0,3].  We  also  have 

Z.V,(  —  )  =  [Z>4(  4-  )4-Zy#(  —  )][Z.S<  4-  )4-Zy,(  —  )]  • 

(A7) 

In  relation  to  branch  0  [Fig.  8(b)]  the  following  re¬ 
lations  also  hold; 
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Z*0(0,I  J  =  +  )+zN<(  -  )J  Zv0(  -  )=  [z*,<+  )+ZW|(  - )] 

Xf e*wZNi(  +  )+ZNi(-)\  X [Z,Vj(  +  )+Z.Vj(  — )]  .  (A9) 

X[e*H'ZN2i  +  )  +  ZNj(-)]  ,  (A8)  Taking  the  ratio  (A8):(A9), 

_ l 

ZsJ.  0, 1  ]  e**[e**^  + 1  )(**%  + 1  )(«*%,  + 1  )ZW>(  -  )Z,V,(  - )  (A10) 

Z^0(  —  )  (  —  ) 


where  yN/  =  [Z*,<  -  )]/[Z^(  4- )]  for  any  i.  From 
(A7)  and  taking  the  limit  N,— ►  «,  Eq.  (A  10)  be* 
comes 

y /2=eA‘-r(€**y  + 1  )V(y  + 1  )■* ,  (All) 

where 

y  *  jim  y*,  • 


containing  the  origin  on  one  branch  of  the  Cayley 
tree  [Fig.  8(b)]  and 

6=j(/-i)"-yjv 

is  the  number  of  ways  of  embedding  a  polymer 
containing  the  origin  on  the  Cayley  tree  [Fig.  8(a)]. 
The  above  expressions  for  a  and  6,  which  are  valid 
for  M  >  2,  will  be  derived  at  the  end  of  this  appen¬ 
dix.  If,  for  convenience,  we  introduce  the  follow- 


Note  that  here  we  have  assumed  the  independence 

ing  new  variables, 

upon  the  particular  site  i.  This  assumption  is 
equivalent  to  neglecting  the  surface  effects. 

(A  17) 

The  monomer  density  ^  in  the  thermodynamic 

P\—Z~'y . 

(A18) 

limit  is  given  by 

(A19) 

4*  lim  Zn{Jr)/Zfj  . 

(A12) 

R  =b/a 

From  Eqs.  (A3)-(A7) 

_ UC-" 

(A20) 

l] 

(y  +  lr  ) 

[(/  —  1 )(/  — 2)  +  2(/—  1  )/M\  ’ 

*~l/| 

(A13) 

Eqs.  (A  15)  and  (A  16)  become 

From  (All) 


♦«- 


{y(efWy  + 1) 


Xy(e*"V  +  D+.H  +  1 


(A  14) 


The  above  calculations  can  be  easily  generalized  to 
chains  made  of  M  monomers  on  a  Cayley  tree  of 
coordination  number  /.  In  this  case  Eqs.  (All) 
and  (A  14)  become 


(fe/q)y  («*">  + 1) 


{6/a>y(e**y  + l)+y  +  l  ’ 

where 

a-i(/-l)*-J[</-2)Af+2] 
is  the  number  of  ways  of  embedding  a  polymer 


(A  1 5) 
(A16) 


(A21) 

*=[/UMi+Mi*>MK**i  +  U+*V*i*  +  1]  •  (A22) 

These  two  coupled  equations  represent  the 
desired  equation  of  state.  Note  that  (A  19)  can  be 
written  in  a  different  way.  In  fact,  from  (41), 

(A  17),  and  (10), 

^  (A23) 

where 

Mp„„(0)--lnfl  2)  +  2]fPu 

(A2A) 

is  the  chemical  potential  of  the  polymers  in  ab¬ 
sence  of  solvent,  WBB  is  defined  in  Eq.  (21b),  and 

K«(0)-(  f/2)WM  (A25) 
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is  the  chemical  potential  of  the  solvent  tn  absence 
of  polymers. 

Now  we  want  to  calculate  explicitly  the  expres¬ 
sion  for  a  =a  ( and  b  =  6(M,/).  According  to 
their  definition  we  have 


a(MJ)  = 


I  /-» 
2 


2  g (M,  )g(Mz) 

W,  +.W ,m,M  -I 


+  (/-l)g(Af  -1)  ,  (A26) 


biM,f)  = 


f 

2 


2  g(MOg(M2) 


+/g<  A/-1).  (A27) 


a  macroscopic  molecule  (infinite  cluster).  Consider 
first  the  case  A#  =2,  p&  =  l.  If  P  is  the  probability 
that  the  ongm  0  [Fig.  8(a)]  occupied  by  a  mono¬ 
mer,  belongs  to  an  infinite  cluster,  then 

p  _  [P(01)  +  P(Q2)  +  P(03)] 

[p(0l  )+p(02)4-/?(03)]  '  (B1) 

PiOi)  is  the  probability  that  there  is  a  dimer  in  (0/) 
and  belongs  to  an  infinite  cluster  (j  =  1,2,3); 
2,P(0i)  is  the  fraction  of  monomers  which  belong 
to  an  infinite  cluster,  p(0i)  is  the  probability  that 
there  is  dimer  in  (Oi)  (i  =  1,2,3),  6  =  ^tp(0i)  is  the 
density  of  monomers.  In  the  limit  of  infinite  sys¬ 
tems  we  have  P(01  )  =  P(02)  =P(03),  and 
p(01)=p(02)=p(03).  Consequently 
P  =P( 01  )/p{ 01 ).  This  quantity  can  also  be  writ¬ 
ten  as 


Here  M ,  and  A/2  >  1  and  g(K)  is  the  number  of 
configurations  of  embedding  a  chain  of  K  bonds  in 
a  given  branch  with  the  origin  fixed  which  is  given 
by 


K>  I 

1,  K  =0 


(A28) 


From  (A26),  (A27),  and  (A28)  for  \f  >  2, 
a(M./)=y(/-l)"-J[</-2)A/+2]  , 

The  ratio 


R=(b /a)  - - J/"11/ _ 

[(/  —  !)(/— 2)  +  2(/-l)/Af]  - 


APPENDIX  B:  PERCOLATION 
OF  INTERACTING  CHAINS 


Here  we  consider  the  connectivity  properties  of  a 
system  of  chains  made  of  M  monomers.  The 
monomers  interact  with  an  effective  interaction 
given  by  Eq.  (10).  Two  chains  are  bonded  if  they 
have  at  least  two  monomers  bonded.  The  probabil¬ 
ity  p$  for  two  nearest-neighbor  monomers  being 
bonded  is  given  by  Eq.  (12). 

We  calculate  now  in  the  Cayley-tree  approxima¬ 
tion  the  gelation  threshold  4f(D,  i.e.,  the  mini¬ 
mum  density  of  monomers  above  which  a  nonzero 
fraction  of  monomers  are  bonded  together  to  form 


P^\-QlQyQAQ,  (B2) 

Qi  is  the  probability  that  there  is  no  infinite  clus¬ 
ter  in  the  branch  (2.6)  — (2.7)  under  the  condition 
that  site  0  is  occupied  by  a  monomer. 

Also  in  the  limit  of  infinite  system  we  have 
Let  us  calculate  Q2  =  Q, 

Q  m  1  occup)+p2(0  °ccup)(?3  ,  (B3) 

where  p2(0  occup)  is  the  probability  that  site  2  is 
occupied  by  a  monomer  under  the  condition  that  0 
is  occupied  by  a  monomer  belonging  to  a  different 
chain, 

p2(o  occup) = ( n0,n2(  +  )>/(n01  > 

=  ( rioi  n  2^  *+•))/[(  rioi  n2(  +  ) ) 

+  <n0in2(-)>] . 

(B4) 

where  n0i  is  the  projector  operator  on  states  in 
which  the  polymer  occupies  position  (01).  n2(  4- ) 
and  n2(  —  )  are  the  projectors  operators  on  states  tn 
which  site  2  is,  respectively,  occupied  or  not  by  a 
monomer.  The  angle  brackets  stand  for  the  ther¬ 
mal  average.  Using  the  same  approach  as  in  Ap¬ 
pendix  A  we  have 

<  n0in2(  +  > )  =  ( 1  /Z.v  )Z.v0[0, 1  ]e^Z.v,l  4- )  . 

(B5) 

( n0,n2(  — )) =( i /Z\ )z.v0[o, i ]Zv,(  — ) ,  (B6) 

where  Zv0(0, 1  ]  is  the  partition  function  of  ti.e 
branch  (0,1)  — (0.3)  [see  Fig.  8(a)],  under  the  condi¬ 
tion  that  there  is  a  dimer  in  (01).  Zv,(  +  )  and 

Z.v2(  —  )  are  the  partition  functions  of  the  other 
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branch  with  origin  in  2  under  the  condition  that 
site  2  is,  respectively,  occupied  or  empty. 

From  lB4)  follows 

p2(0  occup)  =  epwy/(e0wy  -hi) 

(B7) 

Here  y,  ^lt  and  z  have  been  defined  in  Appendix 
A.  In  the  case  of  general  Mt  Eq.  (B2)  becomes 

P  =  \-Q?,  (B8) 

where /=Af(/-2)  +  2.  And  Eq.  (B3)  becomes 

Q=*\-a+aQ?~'  (B9) 


Q  =  \-apB+apBQ?~'  ,  (BIO) 

which  is  the  same  as  Eq.  124)  for  M  =  1  and  h  =0 
with  an  effective  coordination  number  /.  The 
gelation  threshold  will  therefore  be  given  by 

at  =  \/pB{f-\)  ,  (BID 

namely 

AD*  =2/[ ~  1  ]  »  (fii2) 

which  substituted  in  Eq.  (A22)  gives  the  gelation 
threshold 


and  a  Wc  introduce  a  probability  pB 

that  there  is  a  cross  link  between  the  nearest- 
neighbor  monomers.  Then  Eq.  (B9)  becomes 


*t-Jtz2(/-l>/l(n+*M/-l>  —  l]z*  +  (/-  2):j  , 


where/-  !=(/- 1  )pB. 
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